Ternary iron selenide Ko.sFei 6Se2 is an antiferromagnetic semiconductor 
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We have studied electronic and magnetic structures of Ko.8+xFei.6Se2 by performing the 
first-principles electronic structure calculations. The ground state of the Fe- vacancies ordered 
K0.8Fe1.6Se2 is found to be a quasi- two-dimensional blocked checkerboard antiferromagnetic (AFM) 
semiconductor with an energy gap of 594 meV and a large ordering magnetic moment of 3.37 fiB 
for each Fe atom, in excellent agreement with the neutron scattering measurement. The underlying 
mechanism is the chemical-bonding-driven tetramer lattice distortion. Ko.8+xFei.6Se2 with finite x 
is a doped AFM semiconductor with low conducting carrier concentration which is approximately 
proportional to the excess potassium content, consistent qualitatively with the infrared observation. 
Our study reveals the importance of the interplay between antiferromagnetism and superconductiv- 
ity in these materials. This suggests that K0.8Fe1.6Se2, instead of KFe2Se2, should be regarded as 
a parent compound from which the superconductivity emerges upon electron or hole doping. 

PACS numbers: 74.70.Xa, 74.20.Pq, 74.20.Mn 



The recent discovery of high-Tc superconductivity 
in potassium intercalated FeSe[l] and other iron-based 
chalcogenides[2„i3|] has triggered a surge of interest for the 
investigation of unconventional superconducting pairing 
mechanism. In particular, it was found that the super- 
conductivity in these materials coexists with a strong an- 
tiferromagnetic (AFM) order with an unprecedentedly 
large magnetic moment of 3.31 /i^/Fe formed below a 
Neel temperature of 559i^ [1,|5[, meanwhile the conduct- 
ing electron concentration is extremely lowQ. Unlike the 
collinear 0-0 or bi-collinear[l0-12] AFM order observed 
in the parent compounds of other iron-based supercon- 
ductors, the neutron observation found that these mate- 
rials have a blocked checkerboard AFM order 

These compounds have the ThCr2Si2 type crystal 
structure, as shown in Fig. 1(a), isostructural with 122- 
type iron pnictides 0. From the latest X-rays, trans- 
mission electron microscopy, and neutron scattering mea- 
surements, it was suggested that the composition of the 
K-intercalated FeSe superconductors K^Fe2-xSe2 is close 
to KyFei.6Se2 with a fivefold expansion of the parent 
ThCr2Si2 unit cell in the a6 plane, namely a \/b x 
Fe- vacancies superstructure [H, which corresponds 
to one- fifth Fe- vacancies. In an Fe-Fe square lattice, 
an Fe atom nearby a vacancy is 3-Fe-neighbored, or 2- 
Fe- neighbored, or even 1-Fe-neighbored. Topologically 
a square lattice composed only by 3-Fe-neighbored Fe 
atoms exists only in one arranging order but with two 
different kinds of chirality, right-handed or left-handed 
rotations [15], as shown in Fig. 1(b) and (c), respectively. 
This is exactly a x superstructure. 

The peculiar properties observed in K-intercalated 
iron-based chalcogenide superconductors indicate that 
these materials represent a special limit where the inter- 
play between antiferromagnetism and superconductivity 
plays an important role in the formation of supercon- 
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FIG. 1: (Color online) K^Fea.Se2 with the ThCr2Si2-type 
structure: (a) a tetragonal unit cell containing two formula 
units without any vacancy; (b) and (c) are schematic top view 
of an Fe-Fe square layer with one- fifth Fe- vacancies (x = 1.6) 
ordered with the right- and left- chirality, respectively. In this 
>/5 X lattice, each Fe is coordinated with 3 neighboring 
Fe atoms. The squares enclosed by the solid lines denote the 
unit cells. The filled circles denote the Fe atoms while the 
empty circles denote the Fe- vacancies. After optimizing the 
structure by the energy minimization, the Fe atoms go into 
a blocked distribution from a uniform square distribution, in 
which the four closest Fe atoms form a block, represented 
by dashed (blue) circles in (b) and (c), whose bond distance 
is shorter than that of the uniform Fe square lattice (green 
dots). Correspondingly, the angle decreases from 180° to 
176° - 178°. 



ducting pairs. To reveal the hidden physics, we have 
performed the first-principles electronic structure calcu- 
lations for Ko.8+xFei.6Se2 with x = and ±0.1. We find 
that the ground state of K0.8Fe1.6Se2 is an AFM semi- 
conductor with a blocked checkerboard long-range AFM 
order (Fig. 3(c)) and an energy gap of 594 meV, and 
Ko.8±o.iFei.6Se2 are doped electron- or hole- type AFM 
semiconductors. K0.8Fe1.6Se2 is thus an essential par- 
ent compound so that the observed superconductivity in 
K^Fe2-xSe2 is based on the doped AFM semiconductors. 

In our calculations the plane wave basis method was 
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FIG. 2: (Color online) Electronic structure of K0.8Fe1.6Se2 
with one-fifth Fe- vacancies ordered in a way as shown by Fig. 
1(b) or (c) in the nonmagnetic state: (a) the band structure; 
(b-i) are the Fermi surface sheets due to the bands crossing 
the Fermi energy denoted by (1-8) in (a), respectively. The 
corresponding Brillouin zone is shown in Fig. 4(b). The Fermi 
energy sets to zero. 



used. We adopted the generalized gradient approxima- 
tion (GGA) with Perdew-Burke-Ernzerhof formula [16] 
for the exchange-correlation potentials. The ultrasoft 
pseudopotentials [17] were used to model the electron- ion 
interactions. The kinetic energy cut-off and the charge 
density cut-off of the plane wave basis were chosen to be 
800 eV and 6400 eV, respectively. A mesh ofl8xl8x8 
k-points were sampled for the Brillouin-zone integration 
while the Gaussian broadening technique was used in the 
case of metallic states. In the calculations, the lattice pa- 
rameters with the internal atomic coordinates were opti- 
mized by the energy minimization. The optimized tetrag- 
onal lattice parameters in the magnetic states agree ex- 
cellently with the experimental data [5]. 

In the calculations, we adopted a a/5 x ^/Sx 1 tetragonal 
supercell, in which there are two FeSe layers with total 
16 Fe atoms and 4 Fe- vacancies, 20 Se atoms, and 8 K 
atoms and 2 K- vacancies, to represent K0.8Fe1.6Se2. For 
K atoms, there are two inequivalent positions in the su- 
percell. Accordingly, there are a number of inequivalent 
arrangements for the two K- vacancies. From the calcu- 
lations, we find that both the electronic band structures 
and the Fe moments are hardly affected by these inequiv- 
alent arrangements. And the energy difference among 
these different K arrangements is within 1 meV/Fe in 
both spin-polarized and spin-unpolarized cases. This is 
consistent with the experimental observation that the va- 
cancies on K sites are in positional disorder (14)] . 

Figure 2 shows the non-magnetic band structure and 
Fermi surface of K0.8Fe1.6Se2 with one-fifth Fe- vacancies 
ordered in a way as shown by Fig. 1(b) or (c). There are 
8 bands crossing the Fermi level. An important feature 



revealed by the calculation is that due to a substantial 
lattice distortion driven by Fe vacancies, the bond dis- 
tances among four closest Fe atoms around a Se atom 
are reduced. These four Fe atoms will form a compact 
square (we will call it a block below), as shown by the 
light blue circles in Fig. 1(b). The angle ^, which is 
shown in Fig. 1(b), decreases from 180° to about 176°. 
As a result, the Fe-Fe and Fe-Se chemical bonds within 
each block become stronger. This gains an energy of 
about 160 meV/Fe. Such a tetramer lattice distortion 
is not due to the spin-phonon interaction since it hap- 
pens in both non-magnetic and magnetic states. It is 
in fact driven by the chemical bonding, namely this is a 
chemical-bonding-driven lattice distortion. On the con- 
trary, all the structural or lattice distortions previously 
found in the iron-pnictides and chalcogenides are driven 
by the spin-phonon interactions [8], namely magnetism- 
driven lattice distortions. 

Figure 3 shows eight possibly energetically favored 
magnetic orders arranged in a right-chirality lattice (Fig. 
1(b)). From the calculations, if the energy of the non- 
magnetic state is set to zero, we find that the respec- 
tive energies of these eight states, i.e. the ferromagnetic, 
checkerboard AFM, blocked checkerboard AFM, quater- 
nary colhnear AFM, collinear AFM, bi-collinear AFM, 
cell-in collinear AFM, and cell-in bi-collinear AFM states, 
are (-0.070, -0.189, -0.342, -0.232, -0.293, -0.229, -0.254, 
-0.267) eV/Fe for K0.8Fe1.6Se2. The corresponding mag- 
netic moments are (3.29, 2.76, 3.37, 3.10, 3.32, 3.22, 3.20, 
3.26) /i^/Fe. The ground state of K0.8Fe1.6Se2 is thus 
found to have a blocked checkerboard AFM order (Fig. 
3(c)) with a large magnetic moment of 3.37/iB/Fe, in 
agreement excellently with the latest neutron scattering 
measurement [5] . The similar magnetic order was also re- 
ported in a calculation to show metallic KFei.6Se2 0. 

Here we emphasize that it is the chemical-bonding 
driven tetramer lattice distortion mentioned above that 
makes the blocked checkerboard AFM order be the low- 
est in energy. Without this tetramer lattice distortion, 
the energy of the collinear AFM state (Fig. 3(e)) is in 
fact 43 meV/Fe lower than that of the blocked checker- 
board AFM state. The collinear AFM order is the ground 
state magnetic order of the iron-pnictides, driven by 
the As-bridged AFM superexchange interaction between 
two next-nearest Fe atoms, which dominates over other 
exchange interactions. However, in K0.8Fe1.6Se2, the 
chemical-bonding driven tetramer lattice distortion re- 
duces the bond distance within each block so that the fer- 
romagnetic exchange interaction between the two nearest 
Fe atoms inside a block and the AFM exchange interac- 
tions between neighboring two blocks are substantially 
enhanced. It turns out that the blocked checkerboard 
AFM order dominates over all the other magnetic orders. 

Figure 4(a) shows the electronic band structure of 
K0.8Fe1.6Se2 in the blocked checkerboard AFM ground 
state. To our surprise, we find that the compound 
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FIG. 3: (Color online) Schematic top view of eight possible 
magnetic orders in the Fe-Fe square layer with one-fifth Fe- 
vacancies. Among these magnetic orders, the blocked checker- 
board antiferromagnetic order in (c) is the ground state, in 
which the magnetic moments within each block are in paral- 
lel, but the magnetic moments for the neighboring two blocks 
are in anti-parallel. The Fe spins are shown by arrows. The 
empty circles mark the Fe vacancies. The squares or rectan- 
gles enclosed by the dashed lines denote the magnetic unit 
cells. 



K0.8Fe1.6Se2 is an AFM semiconductor with an energy 
band gap as large as ~ 594 meV. This is different from 
the parent compounds of other iron-based superconduc- 
tors, they are either magnetic or non- magnetic semimet- 
als. K0.8Fe1.6Se2 with the checkerboard AFM, collinear 
AFM, or quaternary collinear AFM order is also found 
to be semiconducting. The corresponding energy band 
gaps are 22, 73, and 185 meV, respectively. With any 
other magnetic order, K0.8Fe1.6Se2 is in a metallic state. 

Figure 5(a) shows the total and orbital-resolved par- 
tial density of states (spin-up part) for K0.8Fe1.6Se2 in 
the blocked checkerboard AFM ground state. We find 
that the states around the top valence bands are com- 
posed of both Fe-3d and Se-4p orbitals with the same 
peak energies. This is consistent with the localized fea- 



0.8 



% 0.4 

^ 0.0 

LU 




-0.4 




(b) 




r X M rz A R z 



FIG. 4: (Color online) (a) Electronic band structure of Fe- 
vacancies ordered K0.8Fe1.6Se2 in the ground state with a 
blocked checkerboard antiferromagnetic order (Fig. 3(c)). (b) 
Brillouin zone. Here the top of the valence band is set to zero. 
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FIG. 5: (Color online) K0.8Fe1.6Se2 in the blocked checker- 
board antiferromagnetic ground state: (a) total and orbital- 
resolved partial density of states (spin- up part); (b) projected 
density of states onto the five Fe-3d orbitals at an Fe atom. 
Here the top of the valence band sets to zero. And x-axis is 
along the two next- nearest neighbor Fe-Fe. 



tures of these bands, enhanced by Fe-Fe and Fe-Se chem- 
ical bonding due to the tetramer lattice distortion. On 
the contrary, the lower conduction bands are mainly con- 
tributed by the Fe-3d orbitals. 

By projecting the density of states onto the five 3d or- 
bitals of an Fe atom in the blocked checkerboard AFM 
ground state (Fig. 5(b)), we find that the five up-spin 
orbitals are almost fully filled. This suggests that the 
crystal field splitting induced by Se atoms is small, sim- 
ilar to the case of iron-pnictides [8] . The large magnetic 
moment formed around each Fe results from the Hund's 
rule coupling. The down spin orbitals are partially filled 
by dz2, dxy, and dyz orbitals. Such anisotropy among the 
five down spin orbitals in K0.8Fe1.6Se2 is enhanced by the 
ordered Fe- vacancies. 

We have verified the above results using the virtual 
crystal approximation (VGA) for K atoms. We then re- 
placed K0.8 by K0.7 and K0.9 respectively by using VGA 
to simulate the hole-doping and electron-doping effect 
upon the AFM semiconductor K0.8Fe1.6Se2. As shown 
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in Fig. 6, the electronic band structure of hole or elec- 
tron doped Ko.8±o.iFei.6Se2 is almost the same as the 
one for the un-doped case (Fig. 4). The magnetic mo- 
ments of Fe ions are also hardly changed by doping. 
The doped electrons and holes occupy the bottom con- 
duction bands and the top valence bands, respectively. 
The volumes enclosed by the Fermi surfaces are found 
to be 0.529 holes/cell and 0.642 electrons/cell, namely 
5.00 X 10^°holes/cm^ and 6.06 x 10^°electrons/cm^, re- 
spectively. In principle, every excess Kq.i can supply 1 
electron/cell to the material. Thus both electron and 
hole doping upon K0.8Fe1.6Se2 are very efficient, simi- 
lar to conventional semiconductors. But the total carrier 
concentration is very low, in comparison with other iron- 
based superconductors. This agrees with the infrared 
measurement [i^ . 

We have calculated the formation energy of 
Ko.8+xFei.6_a,/2Se2, defined by (0.8 + x)Ek + (1.6 - 
x/2)EFe + 2^56 — where Ek, Epe, and Ese are the 
atomic energies of K, Fe, and Se, respectively, and E^ 
is the total energy of Ko.8+xFei.6-x/2Se2. For x = 0, 
0.2 and —0.4, corresponding to K0.8Fe1.6Se2, KFei.5Se2, 
and K0.4Fe1.8Se2, we find that the formation energies 
are 16.652, 16.564, and 16.632 eV, respectively. Here 
these stoichiometric compounds are in a balance of 
chemical valence, in which K, Fe and Se are in 1+, 2+, 
and 2- valence states, respectively. The ground states 
of Fe-vacancies ordered K0.4Fe1.8Se2 and KFei.5Se2 are 
in a collinear AFM order [9]. The formation energy 
data indicate that K0.8Fe1.6Se2 is energetically more 
favorable when K atoms are intercalated into FeSe. Thus 
K0.8Fe1.6Se2 is more likely to be the parent compound 
of potassium intercalated FeSe superconductors upon 
doping. The superstructure \/5 x \/b with Fei.6Se2 
layers may thus be the framework or backbone of 
superconductors Ky¥e2-x^^2- 

We have also calculated the electronic structures of 
A0.8Fe1.6Se2 (A=Rb, Cs, or Tl) and found that they are 
all AFM semiconductors. The energy band gaps of these 
compounds are 571, 548, and 440 meV, respectively. 

In conclusion, we have performed the first principles 
calculations for the electronic structure and magnetic 
order of KyFe2-x^^2' We find that the ground state 
of K0.8Fe1.6Se2 is a quasi- two-dimensional antiferromag- 
netic semiconductor with an energy gap of 594 meV. 
The ground state magnetic order is a blocked checker- 
board antiferromagnetic order with a large magnetic mo- 
ment of 3.37/iB, in agreement with the experimental mea- 
surements. The underlying mechanism is the chemical- 
bonding-driven tetramer lattice distortion. Doping elec- 
trons or holes into K0.8Fe1.6Se2 by intercalating more or 
less K to the compound almost do not change the band 
structure as well as the blocked checkerboard antiferro- 
magnetic order. Our result suggests that K0.8Fe1.6Se2 
serves as a parent compound which becomes supercon- 
ducting upon electron or hole doping. 
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FIG. 6: (Color online) Electronic band structures and 
Fermi surfaces in the blocked checkerboard antiferromagnetic 
ground state with the left chirality: (a-c) for K0.7Fe1.6Se2 
and (d-f) for K0.9Fe1.6Se2, to represent the hole-doping or 
electron-doping upon K0.8Fe1.6Se2, respectively. The Fermi 
energy sets to zero. 
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